Characterizing electron entanglement in multiterminal mesoscopic conductors 
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We show that current correlations at the exit ports of a beam splitter can be used to detect elec- 
tronic entanglement for a fairly general input state. This includes the situation where electron pairs 
can enter the beam splitter from the same port or be separated due to backscattering. The proposed 
scheme allows to discriminate between occupation- number and degree-of- freedom entanglement. 

PACS numbers: 03.65.Ud,03.67.Mn,73.23.-b 



c3 



-a 
c 

o 
o 



> 

On 

in 

(N 

o 

o 

-i— > 

c3 



i 

C 

O 
o 



X 



Generation, manipulation and detection of entangled 
electrons is central for realizing integrated solid-state 
quantum computers. Among the several possibilities, a 
lot of attention has been devoted to the study of en- 
tanglement in multiterminal mesoscopic conductors (see 
Refs. [l|, [2| for a review). In this context, it was shown 
that entanglement between spatially separated electrons 
can be detected by means of a beam splitter (BS) Q. 
Indeed the BS, allowing the incoming (and possibly en- 
tangled) electrons to be interchanged, gives rise to two- 
particle interference effects. As a result, the symmetry 
of the incoming state influences the current-noise corre- 
lations at the exit ports. Bunching (enhanced) and anti- 
bunching (suppressed) behavior in the shot noise were 
predicted for spin singlet and triplet entangled states, 
respectively Q. The role of entanglement was later an- 
alyzed in the whole probability distribution of the cur- 
rent fluctuations (being the noise power only its second 
moment) [4]. Further analysis were subsequently per- 
formed in the presence of spin-orbit coupling and for 
states generated in an Andreev double-dot entangler [□]. 
More recently, Burkard and Loss [7[ found a bound for 
the entanglement of arbitrary mixed spin states through 
shot-noise measurements by applying a reduction map- 
ping into Werner states Q. We generalized this result 
to multi-mode input states by introducing an electronic 
Hong-Ou-Mandel interferometer [9(. 

All the previous analysis rely on the assumption that 
only one electron per port enters the analyzer. Most of 
the electronic entangler devices proposed so far P, 
however, generate states having a finite probability am- 
plitude that two electrons enter the analyzer at the same 
input port This gives rise to two distinct forms of 

entanglement: occupation-number and electronic degree- 
of-freedom entanglement HI 12|. Under this generalized 
initial condition the analysis of the entanglement is com- 
plicated by super-selection rules (SSR) induced by parti- 
cle number conservation [H, [ll , 12 , 13| . 

In this paper we present a detection strategy which 
addresses this more general situation. As in Ref. [||, it 
is based on the study of current correlations at the exit 
ports of a BS as a function of controllable phase shifts. 
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FIG. 1: Sketch of the analyzer. An external entangler pre- 
pares electron pairs in the input state of Eq. {TJ and injects 
them into ports 1 and 2. The electrons are allowed to enter 
the analyzer from the same port (cases |$ii) and ($22)), from 
different ports (case $12)), or in any superposition of the pre- 
vious cases. Electrons propagating along lead 2 undergo an 
additional (orbital/spin-dependent) controllable phase shift 
(white circle in the figure) before impinging on a 50% beam 
splitter BS. Current correlations are measured at 3 and 4. 



We shall show that, for the whole class of two-particle 
input states, simple data processing of the measured cur- 
rent cross-correlators can be used to address separately 
the various entanglement components. Moreover, we 
also account for the case in which less than two elec- 
trons can enter the interferometer due to backscattering. 
Our scheme of detection is particularly suitable when 
the entanglement is generated through superconducting- 
normal metal structures or multichannel non-interacting 
structures [3, El]. 

The proposed setup is described in Fig. Q] We start by 
considering a pair of electrons prepared by an external 
device (the "entangler" we want to monitor) at a given 
energy E above the Fermi sea and injected into the in- 
terferometer at the input ports 1 and 2. A generic input 
state has the form 

I*) = sin6»(cos0 |$ u ) +sin0 |$ 22 )) + cos6> |$i 2 ) , (1) 

where 6,<p £ [0,7r/2] and \$ij) describe two electrons of 
energy E entering the BS from the z-th and j-th port 
respectively = 1,2). The latter can be written as 
!*«} = £«,/) a{ a {E)a) p {E)\&), with a^jE) being 
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the creation operator of an incoming electron in lead j 
with orbital/spin degree-of- freedom index a and energy 
E, while 



where 17, 1 



and $' JJ 2 = 

,(12) 



q satisfy the normal- 



ization condition Y, a ,(3 l^ojll 2 = 1 and Ea,/3 l^a,/3l 

1/2. In this notation |0) = |0)i (g) |0) 2 is the vac- 
uum representing the filled Fermi sea of the two leads. 
The analysis of the entanglement contained in the state 
|\&), can be done naturally by bipartitioning the system 
with respect to the port labels 1 and 2. For 6 = 0, 
the state (JXJ) reduces to the previously studied situation 
[E 0, @1- There are however many cases (e.g. the An- 
dreev entangler 

BUS) 

where the incoming state instead 
has the form |T]) with 9^0. Our approach allows to an- 
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alyze this situation as well. Following Refs. [11|, [12| , the 
task is to identify both occupation-number and electronic 
degree-of-freedom entanglement. The first one coincides 
with the "variance-EPR" entanglement of Ref. 1111 and 
with the "fluffy bunny" entanglement of Ref. [l2|: It 
is present whenever we have a non trivial superposition 
among terms where each lead supports a different number 
of incoming electrons, i.e. |$n), | ^22) , and |$i2). The 
second one, instead, coincides with the "entanglement- 
EPR" of Ref. [ll[ and originates from the component 
|$i2) of Eq. |T]): It is present when the leads 1 and 
2 possess one electron each and are entangled through 
the electronic orbital/spin modes a. Under the con- 
straints imposed by the particle number SSR it can be 
shown 
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|13| that occupation-number and electronic 
degree-of-freedom entanglement are both quantum in- 
formation resources which are crucial to solve specific 
tasks (examples are classical and quantum data hiding 
and standard quantum teleportation, respectively). 

As shown in Fig. [TJ electrons propagating along lead 
2 undergo an additional (channel/spin-dependent) con- 
trollable phase shift before impinging on a 50% BS [lj|. 
For simplicity we assume that the BS is symmetric and 
does not suffer from backscattering. The BS is described 
by a scattering matrix s^ a ' defined through the relation 

hA E ) = Ej'=i 4?' a J'AE), with j = 3,4 and b j>a {E) 
the annihilation operator of an outgoing electron in lead 
j and channel a at energy E. The scattering matrix is 

defined as = y/T/2 I J ^ lipo 

we have imposed that no channel mixing occurs in the 
scattering process at the BS. In the case of spin entan- 
glement, this can be easily implemented by introducing a 
BS conserving spin. For orbital entanglement, a spatial 
separation of the orbital channels may be necessary. Fi- 
nally, the tunable phases (p a are introduced at the input 
port 2 by some externally controlled parameters 

The (dimensionless) current cross-correlator at the 
output ports 3 and 4 is 



For definitcness 
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is the current operator of the j-th port (j = 3,4). In these 
expressions the average (■ ■ • ) is taken over the incoming 
electronic state, T is the measurement time, and v is the 
density of states of the leads (a discrete spectrum has 
been considered to ensure a proper regularization of the 
current correlations) . C34 is a measurable quantity which 
can be written in terms of the shot noise S34 and the av- 
erage current 13,4 = J dt(I^^(t)) /T through the relation 
C34 = [S34 + hv 13*4]. Differently from S34, the C34 
is linear in the input state of the system [lj|. Therefore, 
for any given density matrix p = ^iPel^e) (*&e\ with 
as in Eq. fl} we obtain C^{p) = J2e P^zi{^t)i where 
the correlators C^i^i) take the form 

C 34 (tf) = [1 + w cos 2 9 + v sin 2 0sin(2</>)]/4 . (4) 

Here, w and v are real quantities satisfying \w\, \v\ < 1. 
They depend upon the interferometer phases ip a and the 
input state parameters as 
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Consider first the case where there is just one incoming 
electron per port, i.e., 9 = and |W) = | <I? 12)- In this 
situation one has = e/{hv) and Eq. gives the 
same Fano factor F34 = S34/ (2eV«3*4) = [w — l]/4 of 
Ref. where it was shown that all separable mixtures 
of states ($12) must exhibit non-negative values of w, i.e. 
non-negative values of Czi{p) — 1/4. 

In the case where 9 is generic, the same result can be 
obtained by relating the correlator C34 to the entangle- 
ment of formation Ef(p) [2(| of the input state. This 
is done by applying generalized twirling transformations 
to map all density matrices for ports 1 and 2 into gen- 



eralized Werner states [2lJ. The latter are constructed 



as in Ref. fo[ by introducing a joint orthonormal basis 
for ports 1 and 2 formed by the states \xk)i ® \Xk)2 an d 
l*Lt'> = (IXfe>i®IXfc'>2±IXfe')i«)|Xfe)2)/\/2withfc < k'. 
Here, the label fc enumerates all configurations with two 
or less particles in each port described, respectively, by 
vectors |0)i, o| ja (£)|0)i, a\ ta (E)a\ iP {E)\0)i and |0) 2 , 

e-^4jE)\0) 2 , e~^+^+^aljE)alp(E)\0hM- 
This, indeed, generalizes the approach of Refs. [3, |S(, 
since we have introduced a larger basis that accounts for 
ports occupied by a different number of electrons. Fol- 
lowing the derivation of Refs. @, Q, it is easy to show 
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that Ef(p) can be lower bounded by the quantity 

kk 1 

through the inequality 

E f {p) > £{W{p)) , 



(6) 



(7) 



with 8{x) = H [\ + y/x{\ - xjj for x <E [1/2, 1] and null 

otherwise [here H(x) — — a;log 2 x — (1 — x) log 2 (l — x)]. 
In our case, Eq. © yields the identity 



W(p) = 1 - 2C 34 (p) . 



(8) 



It follows that C 34 ( j o) < 1/4 implies W(p) > 1/2 and 
hence, from Eq. ([7]), Ef(p) > 0. Therefore, we can con- 
clude that also in the general case in which the two elec- 
trons can enter the same port, values of C34 smaller than 
1 /4 are direct evidence of entanglement in the input state. 

The sign of C34 — 1/4, however, is not sufficient to 
discriminate between occupation-number and electronic 
degree-of-freedom entanglement. In order to do so we use 
the symmetries w({ip a + tt/2}) — w({ip a }) and v({(p a + 

tt/2}) = -«({?«}) to define c£\p) = [Cufa {<p a }) ± 
C 34 (p; {ip a + 7r/2})]/2. Equation (g]) shows that C- 
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34 



depends only on the \<&\2) component of Eq. (TTJ, while 
C34 depends only upon the superposition among |$n) 
and |$ 22 ), i.e., 



a 



^(tf) = [l+w cos 2 6>]/4, 
C^4 } (*) = v sin 2 6»sin(20)/4 



(9) 
(10) 



From the discussion of the 6 = case it follows that 
the presence of electronic degree-of-freedom entangle- 
ment can be detected by finding values of <p a such that 
C$ < 1/4 & Vice -versa, we observe that any value 
of C34 ' different from zero is indicative of occupation- 
number entanglement in the system. Indeed, C34 =/= is 
possible only for 9^0 and <j) ^ 0, 7r/2 which correspond 
to have a non-trivial superposition between terms with 
different local occupation number. Moreover, one can re- 
late C34' to the superselection induced variance (SIV), 
V(V) = 4((*|Af 1 2 |*) - (^liVxl*) 2 ) [ll|. This is a mea- 
sure of occupation number entanglement, where N\ = 
^ a\ a a,i,a is the total number operator of particles in 

port 1. It is easy to show that V(~&) = (4C^' /v) 2 + 
4 sin 2 #(sin 2 9 cos 2 <f> + cos 2 9), which allows to set a lower 
bound for the SIV, namely: V(V) > (4C*^ } ^ 2 



We finally discuss the case in which the entangler suf- 
fers from backscattering, a situation often encountered 
in proposed entanglers. In this case, due to conservation 
of the total number of particles [J HH, the most gen- 
eral input state in ports 1 and 2 is described by a convex 
combination R of density matrices p", p' and |0)(0| hav- 
ing, respectively, two, one or zero impinging electrons. It 



reads R = q" p" + q' p' + (1 - q' - q") |0)(0|, where q" 
is the probability that both electrons enter the ports (no 
backscattering occurred) and q' is the probability that 
only one of them enters the ports. The entanglement of 
formation of the state R can be bounded as before. In 
this case, however, Eqs. © and © give 

W{R) = q"W(p")+q'W(p')^q"W(p") 

= q"[l - 2C 34 (p")] = q" - 2C 34 (R) , (11) 

where we used the fact that W r (|0)(0|) = and the 
linearity of C34, yielding C 3 i{R) = q" C 3 i{p") (notice 
that since p' and |0) have less than two electrons they 
do not contribute to the current cross-correlator). From 
the monotonicity of the function £{x) we finally obtain 



E f (R) >£(q" -2C 34 (R)) 



(12) 



Once the transmission probability q" of the setup is 
known, Eq. (jT3J) provides an estimate of the entangle- 
ment of formation of the system. Notice that, differently 
from the case discussed in the first part of the paper, the 
state R is a mixture of terms with different total num- 
ber of particles. In this case the standard entanglement 
of formation Ef(R) measures the amount of entangle- 
ment needed to create the state R without taking into 
account the constraints posed by particle-number con- 
servation (i.e., considering possible decompositions of R 
including coherent superpositions of states with different 
particle number) . For this reason Schuch, Verstraete and 
Cirac [llj have proposed to replace the quantity Ef(R) 
with the SSR- entanglement of formation. This is defined 
as Ej SSR \r) = mm Pil ip i Y2iPiEf('ipi), where the mini- 
mization is performed over all decompostion of R with 
\ipi) being eigenstates of the particle-number operator. 
By construction, Ej SSRS> (R) is always greater or equal to 
Ef(R). In our case it reads 



E ( / SR \R) = q" E f {p") + q> E f (p') > q" E f {p") , (13) 



where we used Ef(\0}(0\) = 0. The rhs term of this 
expression can now be lower bounded by noticing that 
p" represents a density matrix formed by vectors of the 
form (TTJ. Employing Eqs. ([7|) and ((8]) we find 



Ef SR \R)^q" £{l-2C 3i {R)/q" 



(14) 



where we used, again, the relation C34(i?) = q" C^p") 
by linearity of C34. As in the case of Eq. (fT2|) . once the 
probability q" is known, Eq. (|14|) provides an estimate of 
the SSR-cntanglement of formation of the system. 

In the remaining part of the paper we make a con- 
nection between the results obtained so far and some 
proposed electronic entanglers. The situation we con- 
sider applies naturally to the case of superconductor- 
based entanglers. It has been shown Q that Cooper 
pairs can be injected into normal metal ports 1 and 
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2 through a (double) tunnel barrier in the form of a 
wave packet |W) = J2e\^e)- Here, the states \^e) 
have the form given in Eq. fl}, with the difference that 
the two electrons do not possess the same energy. In 
fact, they have opposite sign with respect to the chem- 
ical potential of the superconductor, instead, such that 

=Y, a , ^4A E HA- E ^- When calculat- 
ing the cross-correlator C34 of Eq. @ one notices that, 
actually, the energy labels E and —E can be effectively 
incorporated into the mode labels a and j3 in |3>y). 
This is possible thanks to the limit T — > 00 in Eq. 
which allows only zero-frequency contributions (lu = 0) 
to the current operators (J3J). Two-particle wave pack- 
ets produced by superconductor-based entanglers can be 
then described by Eq. JT]) [up to a normalization factor] 
once the energy index is accounted for by the electronic 
degree-of-freedom one. Moreover, notice that introduc- 
ing phases (f a that can discriminate between different 
energies would eventually lead to the characterization of 
energy entanglement in the injected Cooper pairs. 

Regarding non- interacting electron entanglers, at zero 
temperature they typically produce states of the form 

l*> = IIbI^b) with I*b> as in E q- ©■ The en- 
ergy product corresponds to a discretization with interval 

5E = h/T in an energy window eV, where V is a volt- 
age difference applied between two electron reservoirs [l[ . 
The state describes N = eV/SE incoming electron 
pairs. In the case N > 1, the cross-correlator @ ex- 
hibits extra terms scaling with N(N — 1) which do not 
allow to separate the occupation-number entanglement 
from the electronic degree-of-freedom one - details shall 
be given elsewhere. However, the problem may be solved 
by introducing small measurement times instead, such 
that eV <C SE and N — > 1. It can be shown that a calcu- 
lation of the the cross-correlator @ in these conditions 
leads to equivalent results. 
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